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Abstract 

Oxygen causes an increase in the longitudinal relaxation rate of tissues through its T1 -shortening effect owing to its 
paramagnetic properties. Due to such effects, MRI has been used to study oxygen-related signal intensity changes in various 
body parts including cerebrospinal fluid (CSF) space. Oxygen enhancement of CSF has been mainly studied using MRI 
sequences with relatively longer time resolution such as FLAIR, and T1 value calculation. In this study, fifteen healthy 
volunteers were scanned using fast advanced spin echo MRI sequence with and without inversion recovery pulse in order to 
dynamically track oxygen enhancement of CSF. We also focused on the differences of oxygen enhancement at sulcal and 
ventricular CSF. Our results revealed that CSF signal after administration of oxygen shows rapid signal increase in both sulcal 
CSF and ventricular CSF on both sequences, with statistically significant predominant increase in sulcal CSF compared with 
ventricular CSF. CSF is traditionally thought to mainly form from the choroid plexus in the ventricles and is absorbed at the 
arachnoid villi, however, it is also believed that cerebral arterioles contribute to the production and absorption of CSF, and 
controversy remains in terms of the precise mechanism. Our results demonstrated rapid oxygen enhancement in sulcal CSF, 
which may suggest inhaled oxygen may diffuse into sulcal CSF space rapidly probably due to the abundance of pial 
arterioles on the brain sulci. 
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Introduction 

Oxygen causes an increase in the longitudinal relaxation rate of 
tissues through its Tl -shortening effect owing to its paramagnetic 
properties, no T2-shortening effect of oxygen is seen [1,2]. Due to 
such effects, magnetic resonance imaging (MRI) has been used to 
study oxygen-related signal intensity (SI) changes in various body 
organs, such as the lungs [3], brain [4], spleen, myocardium, 
subcutaneous fat, kidneys, bone marrow, liver and arterial blood 
[2,5,6]. The paramagnetic effect of deoxyhemoglobin has 
frequently been used in the brain to visualize blood-oxygen- 
level-dependent (BOLD) contrast as functional MRI [7,8], and the 
paramagnetic effect of the oxygen molecule itself has also been 
used to quantify the oxygen content of cerebrospinal fluid (CSF), 
and to visualize oxygen enhancement (OE) of CSF on MRI [9,10]. 
Fluid attenuated inversion recovery (FLAIR) imaging has been 
mainly used to visualize OE of CSF, since oxygen administration 
causes signal hyperintensity in CSF of the subarachnoid space on 
FLAIR [11,12]. 

CSF acts as a physical cushion for the brain, and plays an 
important role in its biological waste disposal. This fluid is known 
to be produced from the choroid plexus in the ventricles, transfers 
to the cisterns and is eventually absorbed by the arachnoid villi, 
after exchanging contents with the interstitial fluid of the brain. 



CSF has been widely visualized using MRI techniques such as MR 
cisternography [13,14], phase-contrast MRI [15,16] and MRI 
with inversion pulse technique [17], which have provided clues to 
various pathological processes occurring in the brain. Changes in 
the oxygen content of CSF are reportedly associated with injury to 
brain tissue [18], but the resting state of oxygen content and 
dynamic changes in oxygen content after oxygen inhalation 
remain unclear. Knowledge of oxygen changes in CSF is also 
important from the perspective of partial volume effects during 
imaging analysis, since most cerebral cortices and vessels are 
surrounded by CSF. 

OE on MR (OEMR) imaging of CSF has mostly been studied 
using FLAIR [1,11,19], Tl value calculation with inversion 
recovery (IR)-sequences [2,10], and dynamic tracking has been 
performed with relatively longer time resolution [10,11]. Shorter 
imaging time without IR will lead to better temporal resolution, 
but sequences without IR have not been utilized. Fast advanced 
spin echo (FASE) is a similar sequence with half-Fourier 
acquisition single-shot turbo spin-echo (HASTE), which has been 
previously used to calculate Tl values with IR (IR-HASTE) 
[2,20]. This study compared FASE to FASE with IR (IR-FASE) in 
terms of the ability to dynamically track OE of CSF. We also 
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compared OE of sulcal GSF (CSFs) with that of ventricular CSF 
(CSFv) in each image. 

Materials and Methods 

Subjects 

The approval of the ethics committee of Kyoto University 
(approval number: C491) and written informed consent were 
obtained. Fifteen healthy volunteers (12 men, 3 women; mean age, 
32 ±6 years) were recruited, and written informed consent was 
obtained from all volunteers prior to enrolment. 

MRI parameters 

IR-FASE and FASE images were acquired using a 3-T MRI 
scanner (Toshiba Medical Systems, Otawara, Japan) using a 1 3- 
channel head coil and the following parameters in a single-slice 
axial acquisition at the level of 15 mm superior to the anterior 
commissure - posterior commissure line. IR-FASE: repetition 
time (TR), 9575 ms; echo time (TE), 48 ms; inversion time (TI), 
1915 ms; matrix, 192x192, 1.37x1.37 mm; slice thickness, 
7 mm; field of view (FOV), 263 x263 mm; flip angle, 90°; 
bandwidth, 977 Hz/pixel; Number of averaging is two with 
additional one TR for echo stabilization. FASE: TR, 4500 ms; 
TE, 48 ms; matrix, 192x192, 1.37x1.37 mm; slice thickness, 
7 mm; FOV, 263x263 mm; flip angle, 90°; bandwidth, 977 Hz/ 
pixel; Number of averaging is two with additional one TR for echo 
stabilization. 

Dynamic OEMR 

Using IR-FASE (28.8 s/image, 38 measurements) and FASE 
(13.5 s/image, 80 measurements) to track OE of CSF, images 
were divided into three phases: 1) pre-oxygen administration (Pre- 
0 2 ), where subjects breathed normal room air (21 % 0 2 ) for 5 min; 
2) 100% oxygen administration, where subjects breathed 100% 0 2 
at a flow rate of 15 L/min for 5 min; and 3) post-oxygen 
administration (Post-0 2 ), where subjects breathed normal room 
air (21% 0 2 ) for 8 min. Oxygen was delivered through a non- 
rebreather mask that was firmly attached to cover the mouth and 
nose of the subject. 

Image analysis 

Images from each subject were segmented into CSF and non- 
CSF components using a trainable segmentation plugin of Fiji 
[21]. The segmented image was further processed as follows: a 
CSFs mask image and a CSFv mask image were created. CSFs 
and CSFv masks were then applied to each image (38 images for 
IR-FASE, 80 images for FASE) and total SI values for each image 
were calculated (Fig. 1). 

Normalization was achieved by setting the mean SI of all 
dynamic images of Pre-0 2 , 0 2 and Post-0 2 as 1000. Maximum SI 
(maxSI) of CSF was calculated for each subject from the start of 
0 2 inhalation. Slope of SI (SI slope ) was calculated using a 
differential function for SI curve and maximum SI slope (maxSI slope ) 
of CSF was then calculated for each subject from the start of 0 2 
inhalation. Approximation methods were determined by selecting 
best R 2 value for each approximation: linear curve fitting 
regression analysis was performed for the Pre-0 2 and Post-0 2 
phases and polynomial curve fitting regression analysis for 0 2 
phase. Mean SI with standard error of CSFs and CSFv signal 
values in both IR-FASE and FASE independently, and R 2 values 
were calculated for each phase. 




CSFv 



Figure 1. The mean image from FASE is segmented to create 
CSFs and CSFv masks, which are later applied to each image 
and total SI values for each image are calculated. The same 
process is performed for IR-FASE images. Note that the cavum velum 
interpositum was not included neirther to CSFv nor CSFs since it 
belongs to the cistern and it is not located at surfaces, 
doiil 0.1 371/journal.pone.01 00723.g001 

Statistical analysis 

Paired Student's t-test was conducted to compare maxSI and 
maxSI s iope values for each volunteer in CSFs and CSFv between 
IR-FASE and FASE, with values of p<0.05 considered statistically 
significant. MedCalc version 12.2.1 software (MedCalc Software, 
Mariakerke, Belgium) was used. 

Results 

IR-FASE images 

Pre-0 2 showed a linear correlation with time (R 2 = 0.06 for 
CSFs, and R 2 = 0.61 for CSFv). With oxygen administration, 
signal values of CSF increase in correlation with time, with a better 
polynomial curve fit for CSFs than for CSFv (R 2 = 0.97 for CSFs, 
and R 2 = 0.76 for CSFv). Post-0 2 signal values decrease with time, 
showing a linear curve fit (R 2 = 0.68 for CSFs, and R 2 = 0.02 for 
CSFv, respectively) (Fig. 2a, c). 

FASE images 

Pre-0 2 signals showed a linear correlation with time (R 2 = 0.88 
for CSFs, and R 2 = 0.77 for CSFv). With oxygen administration, 
signal values of CSF rose in correlation with time, with a better 
polynomial curve fit in CSFs than CSFv (R 2 = 0.94 for CSFs, 
R 2 = 0.72 for CSFv). Post-0 2 signal values decrease with time, 
showing a linear curve fit (R 2 = 0.94 for CSFs, and R 2 = 0.95 for 
CSFv) (Fig. 2b, d). 

Subtraction images 

0 2 minus Pre-0 2 shows a positive SI difference, while Post-0 2 
minus 0 2 shows a negative SI difference in both IR-FASE and 
FASE images (Fig. 3). 

IR-FASE vs. FASE 

CSFs. Values of maxSI and maxSI s i ope were significantly 
higher for IR-FASE than for FASE (p = 0.001 and p<0.0001, 
respectively) (Table 1, Fig. 4a-b). 



PLOS ONE | www.plosone.org 



2 



June 2014 | Volume 9 | Issue 6 | e100723 



Dynamic Oxygen-Enhanced MRI 



IR FASE CSFs 



SI 1000 




FASE CSFs 



-Pre-02 

02 

Post-02 



SI 1000 




-Pre-02 

02 

Post-02 



200 400 600 800 1000 1200 

Time (seconds) 



IR FASE CSFv 



200 400 600 800 1000 1200 



SI 1000 




-Pre-02 

02 

Post-02 



d. 



SI 



Time (seconds) 



FASE CSFv 




-Pre-02 

02 

Post-02 



1000 1200 



1000 1200 



Time (seconds) 



Time (seconds) 



Figure 2. IR-FASE for Pre-0 2 , Post-0 2 and 0 2 phases of OE-MRI of CSFs (a) and CSFv (c), show higher SI with oxygen administration 
in CSFs compared to CSFv. FASE for Pre-0 2/ Post-0 2 and 0 2 phases of OE-MRI of CSFs (b) and CSFv (d), showing higher SI with oxygen 
administration in CSFs compared to CSFv. All data were shown with mean ± standard errors. 
doi:1 0.1 371 /journal.pone.OI 00723.g002 
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Figure 3. IR-FASE image (a), calculated IR-FASE images: "0 2 minus Pre-0 2 " (b), "Post-0 2 minus 0 2 " (c). FASE image (d), calculated FASE 
images: "0 2 minus Pre-0 2 " (e), "Post-0 2 minus 0 2 " (f). Calculated images of "0 2 minus Pre-0 2 " show a positive SI difference (b, e), while calculated 
images of "Post-0 2 minus 0 2 " show a negative SI difference (c, f). Intra-ventricular high signals in "Post-0 2 minus 0 2 " images are assumed to come 
from the highly vascular choroid plexus. 
doi:1 0.1 371 /journal.pone.OI 00723.g003 
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Figure 4. Box-and-whisker plots of IR-FASE vs. FASE maxSI and maxSI S | ope of CSFs (a-b) and CSFv (c-d). CSFs: maxSI and maxSI S | 0pe are 
significantly higher for IR-FASE than for FASE (p = 0.001 and p<0.0001 , respectively) (a, b). CSFv: maxSI and maxSI S | 0pe are significantly higher for IR- 
FASE than for FASE (p = 0.034, and p = 0.0001 respectively) (c, d). 
doi:1 0.1 371 /journal.pone.01 00723.g004 



CSFv. Values of maxSI and maxSI slope were significantly 
higher for IR-FASE than for FASE (p = 0.034, and p = 0.0001, 
respectively) (Table 1, Fig. 4c-d). 

CSFs vs. CSFv 

CSFs showed significant higher maxSI and maxSI s i ope than 
CSFv in both IR-FASE (p = 0.004 and p = 0.0003 for maxSI and 
maxSI slope , respectively) (Table 1, Fig. 5 a-b) and FASE (p<0.0001 
and p<0.0001 for maxSI and maxSI s i ope , respectively) (Table 1, 
Fig. 5 c-d). 



Table 1. MaxSI and maxSI S | 0 p e for IR-FASE and FASE. 



Discussion 

This study demonstrated dynamic tracking OE of CSF on both 
IR-FASE and FASE, since both methods showed positive signal 
increases during 0 2 administration and maxSI data supported 
these findings. Rapid OE after 0 2 administration displayed by 
SIsiope was a l so demonstrated and represented by maxSI s i ope . In 
addition, OE of CSFs was visualized better than OE of CSFv on 
both IR-FASE and FASE. OE differences between CSFs and 
CSFv appeared largely compatible with previous studies of FLAIR 
with 5-min imaging sequence [1] and Tl value calculation with 7- 
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Figure 5. Box-and-whisker plots of CSFs vs. CSFv. Values of maxSI and maxSI S | 0pe are significantly higher for CSFs than for CSFv in both IR-FASE 
(maxSI, p = 0.004; maxSI S | 0pe , p = 0.0003) (a, c) and in FASE (maxSI, p<0.0001; maxSI S | 0pe , p<0.0001) (b, d). 
doi:1 0.1 371 /journal. pone.01 00723.g005 



min imaging sequence [10], particularly in terms of the regional 
differences of Tl in CSF spaces such as basilar cisterns, lateral 
ventricles and cortical sulci. IR-FASE showed more OE of CSFs 
than FASE. Our results support previous reports of differences 
between OE of CSFs and OE of CSFv, but with higher temporal 
resolution than previously described [1,10]. 

CSF is traditionally thought to mainly form from the choroid 
plexus in the ventricles and is absorbed at the arachnoid villi 
[18,22]. While it is also believed that cerebral arterioles contribute 
to the production and absorption of CSF [23], controversy 
remains in terms of the precise mechanisms [24], Capillaries in 
direct contact with the CSF form the blood-CSF barrier, where 
many constituents pass from the intra-arterial environment into 
CSF [25-27]. Since a linear relationship exists between arterial 
partial oxygen pressure (Pa0 2 ) and CSF oxygen tension, 
increasing Pa0 2 levels with 100% 0 2 administration will lead to 
increased CSF levels of 0 2 . Oxygen diffuses into the CSF through 
the blood-CSF barrier, and this exchange occurs more in CSFs 
than in CSFv, due to the abundance of pial vessels on the surface 
of the brain compared to intra-ventricular vessels. The larger 
amount of intra-ventricular CSF might also cause more dilution of 
oxygen. Another cause might be the leaky areas of the blood-brain 
barrier near the pituitary gland, which would further facilitate 
oxygen diffusion into CSFs more than into CSFv. All of these 
mechanisms might contribute to the difference in OE between 



CSFs and CSFv [1]. A close relationship exists between CSF and 
arterial flows. Phase-contrast cine MRI has revealed the age- 
dependence of CSF flow increases and correlations with arterial 
flow [28]. Our results partly supported the idea that cerebral 
vessels are engaged in CSF production, since oxygen rapidly 
transfers to CSF from arteries. 

IR sequences on MRI have been used to study OE of CSF in 
most previous reports [1,11,12]. Such OE was noticed on FLAIR 
images of anaesthetized patients and was initially explained by the 
effects of propofol, which has a similar Tl to CSF that would lead 
to incomplete nulling of the CSF signal [19]. However, later 
studies interpreted the hyperintensity as due to the paramagnetic 
effects of oxygen on CSF, in turn causing T 1 -shortening effects 
[11,12,29]. Arterial oxygen saturation by hemoglobin is close to 
100% in healthy individuals. The partial pressure of dissolved 
oxygen in the blood represents a small fraction of the total oxygen 
concentration (less than 0.3%). The concentration of dissolved 
oxygen in the blood increased after oxygen inhalation and 
dissolved oxygen in the blood will diffuse into tissues according 
to the oxygen pressure gradient [1]. FLAIR is the sequence with 
an IR pulse to nullify the signal of CSF, so changes in Tl 
relaxation time of CSF interfere with suppression of the CSF 
signal. The T 1 -shortening effect has been utilized on FLAIR in 
various situations [30] . As HASTE sequence was used to calculate 
T 1 values with IR [2] , IR-FASE provides a similar T 1 -shortening 
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effect induced by oxygen as FLAIR, so FASE without IR might 
also show a weak T 1 -shortening effect. Since FASE shows T2 
contrast, a T2-elongation effect cannot be excluded, although 
previous reports have claimed no apparent change in T2 
relaxation time with oxygen administration [2]. The precise 
mechanisms remain unclear, and this appears to be the first study 
to dynamically track OE of CSF on FASE images. 

OE-MRI might act as a non-invasive method for delineating 
pathological lesions occupying CSF spaces [31], and to better 
understand changes in oxygen status in the case of CSF-filled 
cystic tumors. OE-MRI might have the potential to non-invasively 
visualize cerebral collateral blood supply in cases of carotid 
occlusive disease through the diffusion of oxygen into CSF. 

Several limitations exist, since these methods require a high 
degree of compliance from the imaged subject; firm fixation of the 
head was conducted in this study to reduce minor motions of the 
head. Second, quantification of OE such as T 1 value calculation 
was not conducted in this study. We focused on higher temporal 
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